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Abstract 

The EGRET results for gamma ray intensities in and near the Galactic Plane have been 
analysed in some detail. Attention has been concentrated on energies above 1 GeV and the 
individual intensities in a 4° longitude bin have been determined and compared with the 
large scale mean found from a nine-degree polynomial fit. 

Comparison has been made of the observed standard deviation for the ratio of these 
intensities with that expected from variants of our model. The basic model adopts cosmic 
ray origin from supernova remnants, the particles then diffusing through the Galaxy with 
our usual 'anomalous diffusion'. The variants involve the clustering of SN, a frequency 
distribution for supernova explosion energies, and 'normal', rather than 'anomalous' diffusion. 
| ft is found that for supernovae of unique energy, and our usual anomalous diffusion, 

clustering is necessary, particularly in the Inner Galaxy. An alternative, and preferred, 
situation is to adopt the model with a frequency distribution of supernova energies. The 
results for the Outer Galaxy are such that no clustering is required. 
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1 Introduction 

> : 

^ ■ Although supernova remnants (SNR) are often invoked as the source of cosmic rays (CR), of 

energy below the 'knee' at about 3 PeV, PQE], there are still a number of imponderables. Here, 
we concentrate on two aspects: 

(i) SN energies are not unique, but have a certain frequency distribution involving some 'super- 
supernovae' |2]- We refer to this model as SSNR. 

(ii) SNR in the Galaxy are not distributed at random, but rather in clusters, in both space and 
time and, preferentially, in spiral arms This is NSNR model ( Normal SNR ). 

In our previous work [5] we used the acceleration model of [2], which assumed that 
particles are trapped for up to 8.10 4 y after the SN explosion and that SN occurred randomly 
in time in the Galaxy (at an average rate of 10~ 2 y _1 ). The distribution in space was drawn 
from the usual average radially symmetrical distribution of SN surface density, with a peak at 
a Galactocentric distance of D=4 kpc and falling rapidly with increasing D in the vicinity of 
the sun ( at D=8.5 kpc ). Here, we go on to make allowance for the clustering of the SN - 
and hence the SNR - in various ways, and, separately, consider a frequency distribution of SN 
energies: the SSNR model. 

The quantity used as an indicator of CR intensity in the Galaxy is the gamma ray 
intensity [B] which has been measured as a function of longitude, latitude and gamma ray 
energy. Most of the present work relates to gamma rays above 1 GeV which are considered to 
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be produced predominantly by CR nuclei (mainly protons) of median energy ~ 40 GeV [Jj. In 
zero'th order of accuracy, the mean CR intensity along a particular line of sight (£, b) is simply 
proportional to the measured gamma ray intensity divided by the column density of target gas. 
In the actual analysis we use a more elaborate method but the basis is the same. 

The method is to study the ratio of the observed gamma ray intensity to the large 
scale average and use that as an indicator of the small scale variability of CR intensity from 
place to place. The variability can then be compared with that predicted by models involving 
SN, clustered and unclustered, normal and with different energies, using the same technique of 
comparing the intensity with the average. 

2 The Basic Data 

Figure 1 gives an example of the data ( from 6 ); it refers to gamma rays with energy above 
1 GeV and having |6| < 2°, i.e., it relates to the Galactic Plane. Identified discrete sources 
have been removed and the workers estimate that less than 10% of the remaining flux is due to 
unresolved sources. 
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Figure 1: Profile of the gamma ray intensity vs longitude, for _B 7 > 1 GeV and |6| < 2°. The results are from 
UJ. The standard deviations on the points due to Poissonian fluctuations are typically ±2.3% in the Inner Galaxy 
and ±4.8% in the Outer Galaxy. The experimental data are averaged over the longitude bin Al = 1°. The smooth 
line is the best fit of the profile by the 9-degree polynomial. 

We have used the basic data for the numbers of gamma rays to evaluate the statistical 
uncertainty in the individual intensities. 

The data have been divided into 'Inner' (£ : 270° - 350°; 10° - 90°) and 'Outer' 
{£ : 90° — 180° — 270°), the innermost region (\£\ < 10°) being omitted because of confusion 
there. 

Two latitude ranges are considered: \b\ < 2.5° and |6| : 2.5°/5°. The gamma rays from 
the individual b- and /-ranges are derived from regions of the Galaxy distributed along the line 
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of sight but we can identify a median D-value for each range. Restricting attention to gamma 
ray production in molecular clouds with the scale height of ~ 70 pc, the corresponding distances 
are ~3 kpc and 1 kpc from the sun, for |6| < 2.5° and \b\ = 2.5°/5° respectively. Allowing for 
production in HI, with its greater scale height, the D-values will be somewhat different ( but 
by not much because more of the gas - associated fluctuations arise from the clumpy H2 rather 
than HI ). 

3 The Analysis 

3.1 The Measured Fluctuations 

The 'figure of merit' adopted is the ratio of the observed gamma ray intensity to the smoothed 
line through the points ( 9-degree polynomial, equivalent, approximately, to a smoothing over a 
20° longitudinal bin for the Inner and Outer Galaxy regions, separately ). The shape of the line 
is determined essentially by the kiloparsec-scale spatial variations of the cosmic ray intensity 
and the target gas. Figure 1 gives the smoothed curve for \b\ < 2°. 

The ratios have been determined for various bins of longitude, Ai, specifically, 1°, 2°, 4° 
and 8°. Figure 2 shows the corresponding frequency distributions of the differences between the 
EGRET data and the fit. The numbers of entries, the means and the rms (a— ) values are given 
in Table 1. It is with the latter rms values that we are mainly concerned. 




Figure 2: The distribution of differences between the EGRET intensities and the polynomial fit for various 
Ai?-values. The numbers of entries, the means and the RMS (a—) values are given in Table 1. 

The appropriate a- value to be used is a compromise between the A£ over which intensity 
changes are expected (all values but with a bias towards degree scales) and the need to get away 
from small A£-values, where bin-to-bin correlations are serious, because of the finite range of 
latitudes (\b\ < 2.5°) over which the averaging is made. 

We consider that Al = 4° is most appropriate and this will be used in the model 
calculations. 

Figure 3 gives the values of a (4°) for the observed data ( after correction for noise ) 
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and for various model calculations, to be described shortly. The values are plotted at increasing 
(from left to right) distances from the Galactic Centre of the median gamma ray production 
regions. 

3.2 Variations in the column density of gas 

In all our models the radial distribution of gas has axial symmetry. It means that for a smooth 
non-fluctuating CR distribution the longitudinal distribution of the gamma-ray intensity has also 
to be smooth. However, we know that the ISM is highly non-uniform and this non-uniformity 
will contribute to the fine structure of the gamma-ray intensity profile. 

A number of facts are relevant, as follows: 

(i) The large scale dependence of N(H) on I is automatically included through the radial distri- 
bution. 

(ii) The 9-degree polynomial fit causes changes on scales of about 51 = 20° ( and above ) to be 
allowed for. 

(iii) Much of the dumpiness of the gas in the ISM comes from molecular clouds and we allow 
for correlation between the SNR and the clouds in the analysis. 

(iv) Results of the type shown in Figure 2 show a dependence of a on Al in the range A/ : 1° — 8° 
which is very similar for all the models and the observed values. 

Although the implication of the above is that the effect of the fluctuations in N(H) 
( both atomic and molecular ) is not likely to be dramatic, it is certainly finite. Its magnitude 
will be greatest in the Inner Galaxy at low latitudes, \b\ < 2.5°, where contributions from the 
'far' Inner Galaxy are important. 

Estimate for both the Inner and Outer regions have been made. For N(HI) the esti- 
mated value of the logarithmic standard deviation is 0.027 using the well known column den- 
sities of HI |Sj. For N(H2), data [Oj give a value 0.026 referred to the whole column density 
( N(HI) + N(H2) ), leading to an overall value of a gas = 0.037. Comparison can be made with 
the observed value of a Q b s for this I, b region: 0.080. Subtraction in quadrature leads to a = 
0.071, i.e. a small, but finite reduction. 

Support for our analysis comes from a study of the EGRET observations and analysis 
published prior to the latest observations used by us [HI El- In this work, predictions were 
made for a model in which the column densities of gas were used together with an inferred 
CR intensity distribution which was correlated with the volume density of gas ( derived from 
the column density and assumed rotation curve ) but 'smoothed' spatially at the kpc level. 
The majority of the excursions in predicted intensity were thus due to the gas column density 
variations alone. The fluctuations of the observed intensities ( E~ > lGeV, \b\ < 2°, \l\ : ±90° ) 
about the prediction have a = 0.06 ± 0.015. This value is quite consistent with ours. 

Turning to regions away from the Galactic Plane in the Inner Galaxy ( |6| = 2.5° — 5° ) 
the corresponding value for HI and H2 is a gas = 0.04. In the Outer Galaxy the values are 
a = 0.025 for \b\ < 2.5° and 0.035 for |6| = 2.5°/5°. It will be noticed that in all cases the 
corrections will be small. 

Our contention is that the major cause of the fluctuations in gamma-ray intensity on 
the 4° scale is the stochastic nature of the SNR in space and time and, probably, in energy 
output. The column density of gas undoubtedly fluctuates with I but the observed gamma-ray 
intensity fluctuations are much greater because of the CR intensity variations along the line of 
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sight. 



(a) 



Normal diffusion 



INNER OUTER 
|b|<2.5° |b|=2.5°/5° |b| = 2.5°/5° |b|<2.5° 



C 

SSNR 
OBS 



NC 

_J I L 

5 6 



SSNR -str 



OBS 



C 6 



North 
South 



SSNR 

SSNR 

C £ North 
OBS «, South c q 

OBS « 



NC -s?r 



NC -sir 



7 8 9 1 O 

Ga lactocentric distance, kpc 



NC 

_i I i 

11 12 



1 3 



Anomalous diffusion 



CnJ 



^3 



(b) 



|b|<2.5° 



INNER 



SSNR A 



OBS 
C 6 

NC A 



|b|=2.5°/5° 



SSNR A 



OBS 



OUTER 
|b| = 2.5°/5° |b|<2.5° 
SSNR A 



SSNR A 



C 6 



North 



8 INC 
Sc 



NC m 



OBS ^ Sout 



NC m 



c 6 



OBS 



NC 4 



6 7 8 9 1 O 11 

Ga lactocentric distance, kpc 



1 2 



1 3 



Figure 3: Summary of standard deviations (a-values) of gamma ray intensities about the best-fit: (a) - normal 
diffusion, (b) - anomalous diffusion. Small downward corrections have been applied to the OBS ( observed )- 
values to allow for experimental noise. Model simulations are: NC - non-clustering, C - clustering, SSNR - with a 
frequency distribution of energies. Mean standard deviations are obtained by averaging 4-8 independent samples. 
Errors are not indicated for clarity of the figures, but the distribution of a for different samples is broad due to 
the stochastic nature of SN explosions. Typical errors for anomalous diffusion are 0.002-(NC), 0.0015-(C) and 
0.014-(SSNR) for the Inner Galaxy, 0.004-(NC), 0.004-(C) and 0.02-(SSNR) for the Outer Galaxy. For normal 
diffusion the errors are by 2-5 times less. The Galactocentric distances are the means from which much of the 
gamma ray intensities are derived. 



3.3 The model predictions 

Our 'usual' model [3] comprises CR production by SNR, as described in paragraph 1, and 
propagation with 'anomalous diffusion'. The SNR were assumed to be all identical and inde- 
pendent in time and space ( i.e. NSNR ), within the constraints of the dependence of the mean 
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SNR density on Galactocentric radius. We denote this as NC ( non-clustering ), and the results 
for this situation are shown in Figure 3. This is the NSNR model, but with no clustering. 

The SSNR model has the same mean energy content per SNR as in the previous case 
but a distribution of energies about the mean, with a frequence distribution of log energy having 
a log standard deviation of 1.2, following The results for this type of super-supernovae are 
now indicated as SSNR in Figure 3. 

Turning to the possibility of the clustering of NSN, our clustering model has the fol- 
lowing idealised form. At any one place there are SN in sets of 10, distributed at random, in 
time, over 10 6 y. The results for this model 'C ( clustered ), are also shown in Figure 3. Here, 
the SN are identical. 

Finally, calculations have been made for the same situation, as indicated, for 'normal 
diffusion' (ND) rather than our standard: 'anomalous diffusion'. 

3.4 Comparison of observations with the models 

3.4.1 General remarks 

It is apparent that the observed a values are higher than expected for the NSN model without 
clustering - i.e. non-clustering (NC). For normal diffusion NC(ND) the discrepancy is bigger 
still. Below we examine various factors which influence our model predictions. 

3.4.2 Mixed primary mass composition 

The model predictions shown in Figure 3 relate to primary protons. However, the observed 
gamma rays are produced by all CR nuclei, i.e. by the mixed primary mass composition. 
Primary nuclei are more efficient in producing gamma rays than protons of the same energy per 
particle due to the higher interaction cross-section and higher number of interacting nucleons in 
the collision. At the same time fluctuations of the gamma-ray intensity created by CR nuclei 
are lower than those for protons, therefore our model predictions should be reduced for the 
mixed primary mass composition. We calculated the fluctuations for the mixed composition 
with relative abundances equal to 0.416 for protons, 0.265 for helium, 0.170 for CNO group, 
0.149 for iron group of nuclei jlOj . The result is that we have to multiply our predicted a- values 
by 0.93, 0.86, 0.85, 0.88 respectively for the increasing galactocentric distances indicated in 
Figure 3. The uncertainty of this correction is within the range of 0.005 - 0.019. As a result of 
the correction the discrepancy between observed and expected a-values for the non-clustering 
case (NC) with a mixed primary mass composition is bigger still. 

3.4.3 Mode of diffusion 

It is relevant to examine the question of 'normal or anomalous diffusion ?'. The mode of diffusion 
has obvious implications for the present work. Calculations have been made for all cases, i.e. 
non-clustering, clustering and SSNR models using normal diffusion, instead of anomalous. The 
result is that the fluctuations and cx-values are lower for normal diffusion in all cases ( see 
Figure 3a compared with Figure 3b ), thus the difference between the non-clustering case and 
observations will be even larger. 
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3.4.4 Super-supernovae 



Examination of Figure 3 shows that the cx-values for SSNR with a distribution of explosion 
energies is higher than for clustering (C) and for the observed fluctuations (OBS). It is due to 
the basic assumption about the fluctuations of the explosion energy. Added to the fluctuations 
caused by the stochastic distribution of SN in space and time these fluctuations result in a rise 
of the total fluctuations. Inevitably, it is possible to 'dilute' the SSNR values by assuming that 
only a fraction of the SNR are of such high energy. According to j|| the fraction of such SSNR 
is equal to 0.2 for SNIbc and 0.1 for SNIIn, but we regard it as a quantity 'to be determined'. 

3.4.5 Clustering of supernovae and the presence of SSNR 

Although at first sight the case for some clustering of SNR at all distances is strong there is a 
complication that must be allowed for. This concerns the correlation of SNR with target gas, 
particularly in molecular form ( principally H2 ). This correlation is on linear scales smaller 
than the kpc-scale already allowed for. Our 4° bin of longitude would correspond to a linear 
scale of 280 pc at a distance of 4 kpc ( the typical distance to the production region in the Inner 
Galaxy for \b\ < 2.5° ). Such a distance would embrace likely distances of molecular clouds 
(MC) from some, at least, of the SNR. 

The magnitude of this correlation is not clear but some progress can be made. The list 
of known SNR has only 3% of the entries mentioning adjacent MC but this fraction will 
increase with increasing SNR-MC distance. Inspection of the Galactic maps of MC [SJ in the 
Inner Galaxy shows a mean separation of about 500pc; thus a correlation factor of about 0.5 
would appear to be indicated. For the general case, it appears that the relation between the 
various quantities for the situation where there is partial clustering of the SN ( coefficient / ) 
and a partial correlation between SNR and molecular gas ( coefficient 0.5 ) is: 

°L - ° 2 9 as = (1 " <?){(! - 5)[f(R c a c ) 2 + (1 - f)(R nc a 2 nc ] + 5(R ssnr a ssnr ) 2 } 

+g [(l - 6) Rc<J 2 c +Rnc<J 2 nc + § (R ssnr a ssnr f] (1) 

where a Q b s is the total, observed (corrected) standard deviation, a gas characterizes the fluctua- 
tions of the column density of the gas ( as deived in §3.2 ), a nc , a c and 

Cssnr 3-re the expectations 

for non-clustering, clustering and SSNR models respectively and the mixed primary mass com- 
position, R nc , R c and R ssnr - ratios of the model to the observed gamma-ray intensities for the 
same three models, g is the fraction of gas in molecular form and 5 is the fraction of SSNR 
among the total SNR. The first term on the rhs of equation (1) relates to SNR in the ordinary 
interstellar medium, the second term - to SNR inside MC. The equation satisfies the lower limit 
Cobs = <?nc when / = S = g = 

In the calculations it is assumed that SSNR are so rare that they are not clustered, 
whereas half of the ordinary SNR (NC) are clustered in MC. The values of g are, for the approx- 
imate median D-values relating to Figure 3: 0.5, 0.3, 0.1 and 0.1. In fact, different models give 
different absolute gamma-ray intensities, but they can be easily normalised to the experimental 
data multiplying the explosion energy by a constant. The value of relative fluctuations a in this 
operation remains unchanged. Therefore, we adopted all R nc = R c = R ssnr = 1. 

The fraction / is therefore the only unknown variable and it can be calculated from 
equation (1). The results are shown in Figure 4. 
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Figure 4: The fraction of the clustered SNR as a function of galactocentric distance for normal ( o ) and 
anomalous ( * ) diffusion and for 3 SSNR fractions: (a) - 8 = 0, (b) - 8 = 0.15, (c) - a = 0.2. The lines connecting 
the points are 3-degree polynomial fits to the data. The abscissae for the two curves are slightly displaced with 
respect to each other to avoid an overlap of the errors. For 8 — 0, the first point ( extreme left ) for anomalous 
diffusion is at / = 3.21. 



4 Discussion 

The situations for the two models can be considered in turn. 

For the NSNR model, the model used by us in much of our work, considerable clustering 
is needed in the Inner Galaxy. Indeed, with / > 1 the clustering must exceed that adopted by 
us, viz. 10 SN per 10 6 y coincident in position. Although not impossible, it is unlikely that 
the necessary higher number ( 20-30 ? ) is present, even in the congested region of the Giant 
Molecular Ring at D Akpc. 

Instead, we incline towards the SSNR model. Inspection of Figures 4b and 4c shows 
that in the important Inner Galaxy region, with 5 = 0.2, the bulk of the work is done by SSNR 
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and little by SNR. A value of 5 in the range 0.15 to 0.2 is indicated, with 0.15 being more 
physical. In the Outer Galaxy , the likely range is 0.0 to 0.15. 

Certainly, the values of / shown in Figure 4 have big systematic and statistical errors 
due mainly to the simplified model of clustering and big fluctuations caused by the stochastic 
nature of SN explosions. The estimates of these errors are indicated in Figure 4. 

5 Conclusions 

The fluctuations in the intensity of gamma rays above 1 GeV in and near the Galactic Plane 
give information about the mode of production of cosmic rays in their sources, and their mode 
of propagation. 

We have adopted alternative models of both sources ( supernovae of a variety of energies 
and SN of unique energy but with varying degrees of 'clustering' ) and diffusion ( our 'standard' 
anomalous diffusion, and normal ). For 'normal diffusion', although at D ~ 5.5 and 11.5 kpc 
it is possible to achieve 'reasonable' results, for 7 < D < 11 kpc ( the region nearest the sun ) 
this is not possible. Whatever the value of 5, the fractions are unreasonably high. Perhaps this 
argument can be used to point against 'normal' diffusion being applicable for CR of the energies 
in question ? 

Bearing in mind our general arguments that clustering is more common in the Inner 
Galaxy it would appear that <5 ~ 0.15 is favoured ( Figure 4b ). We would then have a constant 
fraction of SN of the SSN variety, with none of them clustered, but the normal SN clustered 
only in the solar vicinity and the Inner Galaxy. 

It seems that the solution put forward is astrophysically reasonable, for the following 
reasons. Concerning standard SNR, these come from stars of modest mass produced in Molec- 
ular Clouds. These clouds are composed of individual clumps in which star production ( and 
subsequent SN ) occurs. In the Inner Galaxy, where the overall MC masses are bigger |12j there 
will be more clusters than in the Outer. For SSN, the fractional number per clump of the re- 
quired massive progenitor stars is probably the same in the Inner and Outer Galaxy. In view of 
continuing problems with theories of star formation [ T.J.Millar, private communication ] even 
this conclusion cannot be regarded as completely firm. 
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Table 1 







No. of entries 


Mean 


o 


1 11 11 QV 


1° 


i-U-L 


n nrwQ^i 


087 




2° 


80 


0.00225 


0.086 




4° 


40 


0.00175 


0.081 




8° 


20 


0.000500 


0.074 


Outer 


1° 


180 


-0.00272 


0.098 




2° 


90 


-0.00156 


0.096 




4° 


45 


-0.00167 


0.091 




8° 


23 


+0.00022 


0.081 



Table 1. Results for the distribution of the differences between the EGRET intensities and the 9- degree polynomial 
fit: nos. of entries, means and RMS (<?-) values. All values relate to the logarithm of the intensity. Latitude 
range, \b\ < 2.5°, Ej > 1 GcV. The histograms are given in Figure 2. 
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